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At zero temperature the lowest part of the spectrum of the QCD Dirac operator is known to consist 
of delocalized modes that are described by random matrix statistics. In the present paper we show 
that the nature of these eigenmodes changes drastically when the system is driven through the finite 
temperature cross-over. The lowest Dirac modes that are delocalized at low temperature become 
localized on the scale of the inverse temperature. At the same time the spectral statistics changes 
from random matrix to Poisson statistics. We demonstrate this with lattice QCD simulations using 
2 + 1 flavors of light dynamical quarks with physical masses. Drawing an analogy with Anderson 
transitions we also examine the mobility edge separating localized and delocalized modes in the 
spectrum. We show that it scales in the continuum limit and increases sharply with the temperature. 



I. INTRODUCTION 

The spectrum of the QCD quark Dirac operator is a 
quantity not directly accessible by experiments but it 
contains essential physical information concerning the be- 
havior of strongly interacting systems. The most well- 
known example of that is the Banks-Casher relation im- 
plying that a non- vanishing Dirac spectral density at zero 
indicates the spontaneous breaking of chiral symmetry 
Another prominent example is that in the interme- 
diate volume so called "epsilon-regime" the lowest part 
of the Dirac spectrum is described by Random Matrix 
Theory (RMT) which makes it possible to extract the 
low-energy constants of chiral perturbation theory from 
the eigenvalues of the Dirac operator (see e.g. Q and ref- 
erences therein). More recently the Dirac spectrum has 
also been used to determine the mass anomalous dimen- 
sion of QCD-like theories with many fermions 

Since at low-temperature the low-lying Dirac spectrum 
has been so extensively studied it is surprising how lit- 
tle is known about it in the high-temperature regime. 
The only solid piece of information above the finite tem- 
perature cross-over is that since chiral symmetry is re- 
stored there the spectral density of the Dirac operator 
should vanish at zero. In principle Random Matrix The- 
ory also has predictions for the spectral statistics at such 
a "soft edge" . However, unlike at zero temperature, at 
high temperature a priori there is no reason to believe 
that the QCD Dirac spectrum is described by this edge 
RMT statistics. Indeed, attempts to verify this numeri- 
cally did not produce fully convincing results @, H| • 

It turns out that random matrix statistics is only one 
of two possible extremes concerning the eigenvalue statis- 
tics. It corresponds to the case of completely delocalized 
eigenvectors occurring only if typical (gauge field) fluctu- 
ations can easily mix eigenmodes nearby in the spectrum. 
The other extreme possibility is localized eigenmodes 
that cannot be mixed by typical fluctuations. In that case 
the spectrum consists of independent eigenvalues obey- 



ing Poisson statistics. Many examples of both types of 
behavior in large linear systems are known both in the 
mathematics and in the physics literature. In fact the 
Bohigas-Giannoni-Schmit conjecture asserts that quan- 
tum systems whose classical counterparts arc chaotic ex- 
hibit random matrix type spectra whereas integrable sys- 
tems after quantization have Poisson-type spectra . An 
example where the same physical system, depending on 
the circumstances, can exhibit both types behavior is An- 
derson localization (loj . In that case the transition of 
single electron states from delocalized ones described by 
RMT statistics to localized states with Poisson statistics 
is driven by impurities in the crystal lattice. 

Already a long time ago the idea was put forward that 
such a transition might also occur in QCD at finite tem- 
perature pd| . Later on numerical studies concluded that 
it is not the case and the Dirac spectrum is described 
by RMT even above the finite temperature deconfining 
and chiral transition [l2[. This conclusion, however, was 
based on the study of full Dirac spectra. It is known that 
in the case of the Anderson transition for weak disorder 
only the states along the band edge become localized and 
states deep inside the band can still remain delocalized. 
Therefore a statistical analysis of full spectra might not 
reveal a localization transition occurring only along the 
band edge. 

With this additional insight the idea of a localization- 
delocalization transition in QCD at T c was revived some- 
time later [l3j]. Both instanton liquid calculations [l4| 
and lattice QCD simulations were J15| done to provide ev- 
idence that such a transition occurs at T c . However, these 
calculations were performed around the critical tempera- 
ture and therefore it was not possible to see clear Poisson 
statistics in the spectrum. For that lattice simulations 
were needed well above T c where localized modes are 
fully developed. Indeed, further support for localization 
in QCD was obtained by a detailed demonstration that 
the lowest two eigenvalues of the overlap Dirac operator 
in quenched SU(2) gauge theory obey Poisson statistics 
[l6| . Finally in the same system a clear transition was 
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observed in the spectrum of the staggered Dirac operator 
from localized Poisson modes to dclocalized RMT modes 
[l7j . The picture emerging from these studies is that 
above the finite temperature transition the lowest part 
of the Dirac spectrum consists of localized modes obey- 
ing Poisson statistics. Higher up in the spectrum there 
is a cross-over to delocalized modes described by random 
matrix statistics. In the meantime, independently, other 
groups also observed the tendency of low modes to be- 
come localized above T c [l8| , although the connection to 
an Anderson-type transition was not made by them. A 
useful account of the spatial structure of low Dirac modes 
and their localization properties can also be found in fl9| . 

So far all the direct evidence for the transition came 
from quenched SU(2) simulations. In the present paper 
we study the question whether full QCD with physical 
light dynamical quarks also exhibits delocalized Dirac 
modes above T c . The question is non-trivial since 
light dynamical quarks suppress the lowest quark modes 
through the quark determinant in the action. Never- 
theless we find that localization also occurs in full QCD 
with quarks of physical masses. We demonstrate this by 
presenting the results of lattice QCD simulations with 
Nf = 2 + 1 flavors of dynamical staggered quarks. We 
also study how the location of the transition within the 
spectrum depends on the physical temperature. 

The paper is organized as follows. At first in Section 
ITT1 we summarize the technical details of our QCD lat- 
tice simulations. In Section UITl we describe the analogy 
between the QCD delocalization-localization transition 
and the Anderson transition. Here we analyze in detail 
the unfolded level spacing distribution that can be used 
to distinguish between localized and delocalized modes. 
We show that in larger spatial volumes the transition 
becomes sharper and most likely in the thermodynamic 
limit it becomes a genuine phase transition. Also in this 
Section we show how to compute the mobility edge sep- 
arating localized and delocalized states in the spectrum 
and analyze localization in terms of the participation ra- 
tio of eigenmodes. In Section lTVl we study how the transi- 
tion scales in the continuum limit. In particular, we show 
that the localization length of localized states is always 
smaller than the inverse temperature. We also demon- 
strate that the properly renormalized mobility edge scales 
in the continuum limit and investigate its temperature 
dependence in the temperature range 1.7T C < T < 5T C . 
Finally, in Section [V] we summarize our results and indi- 
cate further questions. 



II. SIMULATION DETAILS 

At first we summarize the details of our lattice simula- 
tions. We use the Symanzik improved gauge and and the 
two level stout smeared staggered fermion action of Ref. 
[20j with Nf = 2 + 1 flavors. We take the simulation pa- 
rameters from the work of the Budapest- Wuppcrtal col- 
laboration who determined the lattice spacing from the 
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TABLE I. The bare parameters we used in the lattice simu- 
lations; the inverse gauge coupling the light quark mass 
( m ud)i the strange quark mass (m s ) and the corresponding 
lattice spacing (a). 



kaon decay constant Jk and set the bare quark masses 
by requiring the pion and kaon masses to be equal to 
their physical value [2l|. The bare parameters we used 
and the corresponding lattice spacings are summarized 
in Table [J 

To explore the dependence of the localized-delocalized 
mode transition on the temperature and the lattice spac- 
ing wc performed the simulations at three different lattice 
spacings a = 0.06, 0.082 and 0.125 fm and three different 
temporal lattice extensions N t — 4, 6, 8. The physical 
temperature of the system is set by its temporal exten- 
sion as 



T 



1 



(1) 



In this way the lattice parameters we used correspond to 
the physical temperature range of 1.7T C < T < 5T C . 

To understand the nature of the transition it is crucial 
to consider the spatial volume dependence of spectral and 
wave function statistics. For this reason we also repeated 
some simulations on different spatial volumes. The spa- 
tial linear size of the boxes we used were all in the range 
2 fm < L < 6 fm. The details of the parameters of our 
ensembles arc summarized in Table [TTJ 

Since low Dirac modes can potentially be slow modes 
of the system we also checked for autocorrelations. In a 
long run performed on the finest lattice we found that 
the autocorrelation time for the smallest Dirac modes is 
definitely smaller than ten trajectories. To be on the safe 
side on the finest lattice the configurations that we used 
for spectrum calculations were always separated by thirty 
trajectories. Even on the coarsest lattice configurations 
were separated by ten trajectories. 



III. ANDERSON TRANSITION IN THE DIRAC 
SPECTRUM 

Previously it was seen that in the high temperature 
deconfining phase of SU(2) Yang-Mills theory there is a 
transition in the staggered Dirac spectrum from localized 
low modes to delocalized modes higher up in the spec- 
17] . The question we ask here is whether such a 
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transition also occurs in real QCD with light dynamical 
quarks. This is a non-trivial question since light quarks 
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TABLE II. The parameters of the simulations; the tempera- 
ture, the lattice spacing, the spatial and temporal box size, 
the inverse gauge coupling, the number of configurations and 
the number of Dirac eigenvalues computed on each configu- 
ration. 



suppress low Dirac modes through the fermion determi- 
nant in the functional integral and in the quenched case 
it is exactly the lowest Dirac modes that are localized. 

The hallmark of a transition from localized to delo- 
calized modes in terms of spectral statistics is a change 
from Poisson statistics to random matrix statistics in the 
spectrum. Intuitively speaking, localized modes are such 
because they cannot mix with other modes; typical gauge 
field fluctuations do not mix them. Localized modes close 
in the spectrum are peaked at spatially distant locations 
and they are sensitive only to gauge field fluctuations 
there. As a result the corresponding eigenvalues are sta- 
tistically independent and the level spacings obey Pois- 
son statistics. Dclocalizcd modes, in contrast, are mixed 
by typical gauge field fluctuations. Gauge field fluctu- 
ations change several delocalized modes together which 
introduces delicate correlations in the spectrum and as 
a result the eigenvalue statistics is described by Random 
Matrix Theory (RMT). 

The simplest way to detect a transition in the spec- 
trum from Poisson to RMT statistics is to consider the 
so called unfolded level spacing distribution. Unfolding 
is essentially a local rescaling of the eigenvalues to have 
unit spectral density throughout the spectrum. We did 
the unfolding by ordering all eigenvalues in the given en- 
semble and replacing them with their rank order nor- 
malized by the total number of configurations. On a few 
ensembles we also checked unfolding by using local spline 
approximations to the spectral density but there was no 
discrepancy between the two methods of unfolding. 

Since the unfolded level spacing distribution (ULSD) 
is known analytically for both the Poisson and the RMT 
statistics it can be easily used to distinguish between the 



FIG. 1. The unfolded level spacing distribution in different re- 
gions of the spectrum of ensemble C3. The figures correspond 
to the spectral regions 0.15 < Xa < 0.19 (a), 0.29 < Xa < 0.32 
(b), 0.34 < Xa < 0.35 (c) and 0.375 < Xa < 0.385 (d). 
The dashed line indicates the exponential distribution cor- 
responding to the localized (Poisson) case and the dotted line 
indicates the chiral unitary Wigner surmise expected in the 
delocalized (RMT) case. 



two cases. For Poisson statistics the ULSD is a simple 
exponential, 

■Ppoisson(s) = exp(-s). 

In the RMT case the unfolded level spacing distribution 
depends on the universality class which in the case of 
staggered fermions in the fundamental representation of 
the SU(3) gauge group is the Chiral Unitary Ensemble 
(ChUE) 0- The corresponding ULSD is very precisely 
approximated by the chiral unitary Wigner surmise, 

32 / 4 \ 
^ChUE(s) = ~^s 2 ■ cxp s 2 ). (2) 

7T Z \ TT J 



To demonstrate the transition in the spectrum in Fig. [I] 
we plot the ULSD in different regions of the spectrum of 
ensemble C3; 0.15 < Xa < 0.19 (a), 0.29 < Xa < 0.32 
(b), 0.34 < Xa < 0.35 (c) and 0.375 < Xa < 0.385 (d). 
We also indicate in the same Figure the distributions 
corresponding to the Poisson (localized) and the RMT 
(dclocalizcd) case. The transition from localized modes 
at the edge of the spectrum to delocalized modes in the 
bulk can be clearly seen. This shows that light dynam- 
ical fermions do not change the picture observed in the 
quenched case earlier and the transition also occurs in 
QCD with quarks of physical masses. 
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A. Analogy with the Anderson transition 



The transition in the spectrum from localized to de- 
localized modes is reminiscent of the Anderson metal- 
insulator transition occurring in conducting crystalline 
solids when impurities are introduced In that case 
in the presence of impurities single electron Bloch states 
along the band edge turn into localized states. In three 
dimensions if the impurity concentration is not too high 
states at the band center can still remain delocalized. 
The boundary between localized and delocalized states 
is known as the mobility edge. Increasing the density 
of impurities pushes the mobility edge further towards 
the center of the band until all the states in the band 
become localized. When the mobility edge passes the 
Fermi energy and the Fermi energy gets into the delo- 
calized part of the spectrum the system has a vanishing 
zero-temperature conductivity. The states that can be 
excited are all non-conducting localized states. 

It has been conjectured that the finite temperature 
QCD transition might be similar to the Anderson transi- 
tion [l3| . Further indications to sup port this picture were 
obtained from instanton liquid [14| and lattice QCD sim- 
ulations (l5j . 

We now sketch the analogy between the spectrum of 
the one electron Hamiltonian in disordered media and 
that of the QCD Dirac operator. Due to the symme- 
tries of the QCD Dirac operator its spectrum is symmet- 
ric with respect to the real axis. In the continuum and 
also in the case of the staggered lattice Dirac operator, 
the one we use here, the spectrum is purely imaginary. 
That is if the quark mass is zero, otherwise the quark 
mass provides a trivial real part to all the eigenvalues. 
In the chiral limit (zero quark mass) the spectral density 
at zero is proportional to the chiral condensate, the order 
parameter of spontaneous chiral symmetry breaking [l[. 
At high enough temperature chiral symmetry is restored 
and the spectral density at zero vanishes. In that case 
the spectrum has a so called "soft edge" and there might 
even be a gap around zero in the spectrum 0, Q . This 
"edge" of the spectrum at the low end is analogous to 
the band edge in condensed matter systems. 

In the condensed matter literature numerical studies of 
the Anderson transition usually concentrate on the band 
center and locate the critical disorder when states at the 
center become localized. As we will see, in the case of 
QCD the location of the mobility edge is controlled by 
the temperature and there is no analog of the band cen- 
ter since the spectral density continues to be non-zero 
all the way up to the cut-off scale. Therefore, unlike in 
most of the condensed matter literature, here we do not 
attempt to determine the critical disorder where all delo- 
calized states disappear but rather study how the mobil- 
ity edge changes with the temperature. This approach is 
not completely unknown in the condensed matter litera- 
ture either 1221] . 
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FIG. 2. Variance of the unfolded level spacing distribution 
across the transition in the spectrum for ensembles CI, C2 
and C3 that differ only in their spatial volume. The dashed 
horizontal lines at 1 and 0.178 indicate the expected limits in 
the localized (Poisson) and delocalized (RMT) regime. 



B. Second moment of the unfolded level spacing 
distribution 



As can be seen in Fig. [T] the unfolded level spacing dis- 
tribution changes in the spectrum from the exponential 
to the Wigner surmise in a continuous fashion. There 
does not appear to be a sharp mobility edge, A c , separat- 
ing localized and delocalized states. Even in the case of 
Anderson transitions, however, a sharp transition is ex- 
pected only in the thermodynamic limit when the spatial 
size of the system, N s , goes to infinity. This is completely 
analogous to second order phase transitions where a truly 
divergent correlation length and a sharp phase transition 
can only be observed in infinite systems. In principle any 
physical quantity that changes in a well defined way from 
the localized to the delocalized regime can be used to de- 
fine a transition point in the spectrum. If there is a sharp 
transition in the thermodynamic limit then non-analytic 
behavior should appear in all these quantities at a given 
point, A c , in the spectrum. In what follows we will look 
for a quantity that can be used to define the transition 
point. 

Since the first moment of the unfolded level spacing 
distribution is unity by construction, the simplest quan- 
tity to consider is its second moment or variance. In the 
localized case the level spacings are exponentially dis- 
tributed and the variance is 



<* 2 > - <*> S 



1. 



(3) 



while in the delocalized regime the second moment of the 
distribution of Eq. ([2]) can be analytically determined to 
be <7g = ^ — 1. In Fig.[5]we plot how the variance changes 
in the spectrum throughout the transition region. The 
three data sets correspond to ensembles C1-C3 with three 
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FIG. 3. The same as Fig.[T]but zooming in on the transitions 
region. The dotted, dashed and the continuous curve indicate 
the fits of the form Eq. (gj to the N s = 16, 24 and 32 data 
respectively. 



different spatial volumes but otherwise identical param- 
eters. It is apparent that with increasing volume the 
transition becomes sharper. To define a finite volume 
pscudocritical point in the spectrum we locate the inflec- 
tion point of the curve cr 2 (A). To this end we use the 
three-parameter fitting Ansatz 

a 2 (A)=.4{l-tanh(5(A-C))} + ^--l. (4) 

8 

This form ensures the correct limit for large A and yields 
good fits starting already from Aa = 0.2, 0.25 i.e. already 
from below the transition point. Using this Ansatz the in- 
flection point can be easily seen to be at A c a = C and the 
slope of the curve there is AB. These two parameters are 
largely independent of the starting point of the fit as long 
as it starts at smaller values of A than where the inflection 
point occurs. The location of the inflection point turns 
out to be also independent of the volume yielding pseudo- 
critical points C = \ c a = 0.321(4), 0.322(1), 0.324(2) for 
the spatial sizes N s = 16,24,32. The slope at the inflec- 
tion point is AB = 7.6(5), 10.96(36), 13.8(5) for the three 
different spatial sizes and it scales roughly proportionally 
to the linear spatial size of the box. This indicates that 
there might be a genuine sharp transition in the ther- 
modynamic limit. The sharpening of the transition with 
the volume can be better seen in Fig. [3] where we plot 
again the variance of the unfolded level spacing distribu- 
tion but zoom in on the transition region and also show 
the fitted curves with the above described parameters. 

C. Eigenvector statistics 

Besides the spectrum the spatial profile of the corre- 
sponding eigenvectors also contains important informa- 
tion concerning their localization properties. A quantity 



that is widely used in this context is the inverse partici- 
pation ratio (IPR) defined as [23[ 

p^^hM*)! 4 , (5) 

X 

where ip is an eigenvector normalized as 

X>(aO| 2 = l. (6) 

X 

The qualitative physical meaning of the IPR can be eas- 
ily seen by noting that an cigenmode that spreads uni- 
formly in a four-volume v and is zero everywhere else 
has IPR = 1/v. The IPR -1 thus measures the volume 
occupied by an eigenmode. Alternatively one can use the 
participation ratio, PR = IPR -1 ■ V, where V is the total 
volume of the system. This measures the fraction of the 
total volume occupied by the eigenmode. 

The behavior of the IPR and the PR in the thermo- 
dynamic limit can be used to distinguish between local- 
ized and delocalized modes. By definition a part of the 
spectrum consists of localized modes if their average IPR 
remains finite as the volume goes to infinity. This also 
implies that their average PR vanishes in the thermody- 
namic limit. In contrast delocalized modes have vanish- 
ing IPR which usually implies non-vanishing PR. 

In Fig. [4] we plot how the average participation ra- 
tio of cigenmodes changes throughout the spectrum on 
ensembles C1-C4. The four data sets correspond to dif- 
ferent spatial sizes but otherwise identical parameters. 
It is clearly seen that the average PR for the low eigen- 
modes decreases as the volume is increased; the aver- 
age PR here tends to zero in the thermodynamic limit. 
This suggests that these eigenmodes are localized. They 
fill a vanishing fraction of the total box volume in the 
thermodynamic limit. Higher up in the spectrum the 
average PR becomes a volume independent constant of 
order unity, which means that these eigenmodes are de- 
localized. They fill a non-zero fraction of the total box 
volume in the thermodynamic limit. 

IV. CONTINUUM LIMIT AND 
RENORMALIZATION 

We saw that on all the ensembles that we considered 
there is a critical point ( "mobility edge" ), A c , in the Dirac 
spectrum that separates localized and delocalized eigen- 
modes. If the lowest part of the spectrum becomes lo- 
calized that can have a dramatic effect on long distance 
correlators of quark operators and the masses associated 
with them. This happens because in the spectral de- 
composition of the quark propagator, (D + to) -1 , each 
eigenmode is weighted by the inverse of the correspond- 
ing eigenvalue. Therefore the lowest eigenmodes receive 
the largest weight. On the other hand, eigenmodes local- 
ized to a distance scale I have negligible contribution to 
correlators at distance scales much larger than I. 
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FIG. 4. Average PR of the eigenvectors for lattices with the 
same temperature and lattice spacing but with four different 
aspect ratios. N t = 4 



To assess the physical implications of localization in 
QCD we are thus lead to study two questions; 

1. What is the distance scale I on which the lowest 
cigcnmodcs are localized? 

2. How far up in the spectrum (A c ) are the modes 
localized? 

Both of these questions have to be considered in the con- 
tinuum limit as the lattice spacing a — > 0. 



A. Localization length 

To get a rough estimate of the localization scale we 
can consider the inverse participation ratio (IPR) defined 
by Eq. [5] Since the IPR scales like the inverse four- 
volume of the region where the given mode is spread out 
a good measure of the localization scale is provided by 
the quantity 



I 



(IPR--), 



(7) 



where the average is understood over all the eigenmodes 
in a given region of the spectrum. In principle I varies 
through the spectrum but it does not change too much 
within the region of localized modes. To illustrate that, 
in Fig. [5] we show in two representative cases how I 
changes through the spectrum. As can be seen in the Fig- 
ure for the lowest part of the spectrum the localization 
length is almost constant and is independent of the spa- 
tial volume. In this region there might be a small dip in 
I. The dip is generally more pronounced at lower temper- 
atures and/or on finer lattices but on our ensembles the 
total variation of I in the localized regime never exceeds 
20%. At some point in the spectrum I starts to increase 
sharply and becomes strongly volume dependent. This is 
the beginning of the transition to delocalized states. 
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FIG. 5. The localization length (Eq. [7]) of eigenmodes along 
the spectrum for ensembles Dl and D2 differing only in their 
spatial volumes. Both I and the location in the spectrum, A 
are given in lattice units. Statistical errors are not shown as 
they are smaller than the size of the symbols. 



In what follows we define the localization length of lo- 
calized modes with the following simple procedure. The 
localized eigenmodes all have I's between the bottom of 
the dip and the I of the very lowest eigenmodes. This 
interval for ensembles Dl and D2 is indicated by the two 
dashed horizontal lines in Fig. [5] The central value we 
quote for I is always the center of this band and the un- 
certainty is half the width of the band. Compared to 
that, statistical errors are always negligible. 

Having a well defined measure of the localization 
length for localized modes we can now look at how it 
depends on the lattice spacing. Since, as we will see, I 
also depends on the temperature, we choose to compare 
ensembles C, D and E which are almost exactly at the 
same physical temperature of about 400 MeV. Since I 
does not depend on the spatial volume we omit the num- 
bers from the ensemble labels here. In Fig. [6] we plot the 
localization length as a function of the lattice spacing for 
these three ensembles. 

To guide the eye we also included a linear fit to the 
data. Even if the quality of our data does not allow 
a proper continuum extrapolation it can be safely con- 
cluded that the localization length measured in physical 
units does not increase in the continuum limit and it is 
not larger than a few tenth of a Fermi at this tempera- 
ture. 

It is also instructive to see how the localization length 
compares to the most important length scale in the prob- 
lem, the inverse temperature or in other words the tem- 
poral size of the box. In Fig. [7] we plot the localiza- 
tion length in units of the inverse temperature for all the 
ensembles. Different symbols represent data sets corre- 
sponding to different values of the lattice spacing. There 
is a slight trend of the finer lattices producing more lo- 
calized low modes as can be seen also in Fig. [5] The 
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FIG. 6. The localization length of the localized modes as a 
function of the lattice spacing. The three points correspond to 
the same physical temperature of about 400MeV (ensembles 
C, D and E). 
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FIG. 7. The localization length in units of the inverse temper- 
ature as a function of the physical temperature. The different 
symbols represent data obtained at different values of the lat- 
tice spacing. 



scale the system behaves as if it had a gap of order A c . 
The influence of this effective gap on the physics depends 
on how A c , the critical point in the spectrum, scales in 
the continuum limit. 

In Subsection IIII Bl we showed how to determine the 
critical point using the second moment of the unfolded 
level spacing distribution. This procedure yields the crit- 
ical point X c a in dimensionless lattice units for each tem- 
perature and lattice spacing. In what follows we give a 
proposal for defining the continuum limit of this quan- 
tity. Since the critical point is effectively a gap for quark 
modes capable of propagating to long distances it plays 
a role similar to the quark mass that also introduces a 
gap. For this reason we expect A c to be renormalized in 
exactly the same way as the quark mass and the ratio 
of X c a to the bare light quark mass m u da should have 
a proper continuum limit. Moreover this quantity mea- 
sures the relative size of the effective gap for delocalizcd 
modes and the gap for all modes provided by the quark 
mass. 

Another, perhaps more formal, argument showing that 
A c /m u d has a well-defined continuum limit is as follows. 
On the one hand, the pseudoscalar meson correlator 
(P(x)P(0)) is proportional to the matrix elements of the 
square of the Dirac propagator, [(D^ +m)(D+m)}~ 1 . On 
the other hand, for asymptotically large temporal sepa- 
rations t 



i£(PM)P(0)} = Cp P e" 



(8) 



where V s is the spatial volume and Cpp is related to the 
pion decay constant, f n , as [Hf 



(9) 



The right-hand side of this equation is a well-defined 
physical quantity and thus so is the left-hand side. Since 
Cpp is proportional to the inverse quark propagator 
squared, it is also proportional to the inverse of the Dirac 
eigenvalues squared. Thus 



jA c 2 oc C PP m 2 d 



(10) 



important point is that the lowest modes always appear 
to be localized on or below the scale of the inverse tem- 
perature in the whole range of temperatures studied here. 



B. Renormalization and temperature dependence 

of A c 



We have seen that the lowest eigenmodes spread out 
only to a distance scale below the inverse temperature. 
Therefore these modes cannot contribute to quark prop- 
agation on length scales larger than that. As far as 
hadronic correlators are considered above that distance 



also has a well defined continuum limit. For a more gen- 
eral discussion of similar issues with Wilson fermions see 
also [H and [13. 

Data sets C, D and E where we have data at the same 
temperature for all three lattice spacings indicate that 
the lattice spacing dependence in the quantity is 
comparable to its uncertainty at a given lattice spac- 
ing. Therefore in Fig. [8] we plot the temperature de- 
pendence of this quantity for all lattice spacings in the 
same figure. The different types of symbols, correspond- 
ing to data obtained from lattices of different coarseness, 
all lie on a smooth curve confirming that scaling viola- 
tions in this quantity are small. The lowest order poly- 
nomial fit to all the data that yields an acceptable chi- 
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The continuous line is a quadratic poly- 



nomial fit to all the data. The different symbols correspond 
to data obtained using lattices of different coarseness. 
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FIG. 9. The number of localized (below A c ) eigenmodes per 
cubic Fermi as a function of the temperature. Different sym- 
bols refer to simulations at different values of the lattice spac- 
ing. 



that means practically in Fig. [9] wc plotted the average 
density of eigenmodes below A c as a function of the tem- 
perature. Again we collect data obtained at different val- 
ues of the lattice spacing in the same Figure. Apparently 
the number density of localized modes decreases sharply 
as the transition temperature is approached from above. 
In practical terms that means that in order to have for 
instance about twenty localized modes per configuration 
the aspect ratio of the boxes has to be kept between 4-6 
in the temperature range considered here. 

Finally we would like to point out a potential techni- 
cal difficulty in studying the level spacing distribution. 
It is caused by the pairing of staggered Dirac eigenval- 
ues affecting the lowest part of the spectrum when the 
spatial volume is small. It comes about because in the 
continuum limit staggered fermions describe four fermion 
flavors and the spectrum is expected to become fourfold 
degenerate. Therefore approaching the continuum limit 
this degeneracy starts to be formed and peculiar corre- 
lations will appear among members of the would-be con- 
tinuum quartets. 

It turns out that already at the lattice spacings we 
used here the stout smearing of gauge fields coupled to 
the fermions bring the spectrum close enough to the con- 
tinuum behavior that such extra correlations for the low- 
est eigenvalues can be detected. This is first manifested 
in pairwise attraction between consecutive eigenvalues 
which distorts the Poisson statistics [29]- The splitting 
between members of the doublets, however, does not de- 
pend appreciably on the volume while the separation of 
the doublets (eventually quartets) should be inversely 
proportional to the spectral density and thus to the vol- 
ume. Therefore this distortion of the level spacing statis- 
tics is a finite volume effect (c.f. also [30j]). Moreover it 
only affects the lowest part of the spectrum where the 
spectral density is very small. Since the spatial volumes 
we used are large enough the results we presented here 
are not affected by the staggered degeneracy. On the 
other hand, a more systematic study of the level spacing 
statistics of the very lowest staggered Dirac modes might 
reveal interesting facts about how the expected degener- 
acy is formed in the continuum limit. 



squared (x 2 = 1.3) is second order and extrapolates to 
T =170McV at = 0. This value is consistent with 

m ud 

the known location of the finite temperature cross-over 
in QCD [27], HH below which chiral symmetry is broken 
and no localized modes are expected to be present. This 
provides a further consistency check of our results. 

Another physically interesting quantity is the number 
density of localized modes. This quantity is also practi- 
cally important if one aims at detecting localization from 
spectral statistics. In order to be able to observe Poisson 
statistics in the level spacing distribution the volume of 
the system has to be large enough to accommodate sev- 
eral localized eigenmodes per configuration. To see what 



V. CONCLUSIONS 

In the present paper we argued that in QCD above 
the finite temperature cross-over the lowest eigenmodes 
of the quark Dirac operator are localized. The spatial 
localization length is set by the inverse temperature (see 
Fig. [7|) with eigenmodes becoming more "squeezed" at 
higher temperature. At the same time when the temper- 
ature increases the mobility edge, separating localized 
and delocalized modes, is also pushed higher up in the 
spectrum. 

This high temperature behavior of the low Dirac modes 
has to be contrasted with the situation at low tempera- 
ture below the cross-over. In that case chiral symmetry is 
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spontaneously broken and a finite density of cigenmodes 
extends down all the way to zero. As a result, the statis- 
tics of the lowest Dirac modes is described by random ma- 
trix theory. QCD is thus a remarkably rich theory. The 
low-end of the Dirac spectrum can exhibit both possible 
extremes of spectral statistics; maximally mixed modes 
with RMT statistics below the transition and completely 
independent cigenmodes with Poisson spectral statistics 
at high temperature. 

The localization of the lowest Dirac modes can dramat- 
ically suppress hadronic correlators at high temperatures. 
In the eigenmode expansion of the quark propagator the 
lowest part of the spectrum receives the largest relative 
weight. If, as we saw, these modes are localized they 
cannot propagate quarks to long distances resulting in a 
suppression of long-distance correlators. As can be seen 
in Fig. [8] the mobility edge below which all states are 
localized moves steeply up with increasing temperature. 
Already at 2T C it is two orders of magnitude larger than 
the bare light quark mass. This mechanism might help to 
explain the steep rise of screening masses above T c seen 
in lattice simulations |3lj . 

An interesting question is how exactly the spectral and 
wave function statistics change through the mobility edge 
and whether there is any universality in how these quan- 
tities interpolate between the localized and delocalized 
regime. In particular we would like to check whether the 
scale invariance of the inverse participation ratio distri- 
bution observed in Anderson transitions [32j also occurs 
in the QCD transition. We hope to return to this ques- 
tion in a future publication. As far as the Poisson to 
RMT transition in the spectral statistics is concerned 



there does not seem to be a universal understanding of 
how it happens in general but there are several proposals 
that might provide further insight [33|-[3t|. 

It would also be interesting to know what physical 
mechanism drives the transition. Is it possible to identify 
some physical objects in the gauge field background that 
are responsible for the appearance of localized modes? 
In a previous paper some evidence was found that there 
is a correlation between localized modes and local fluc- 
tuations of the Polyakov loop [lij . On the other hand it 
was also argued there that uncorrclated instantons can- 
not play a significant role in this mechanism as their den- 
sity is too low for that. It is, however, still possible that 
instanton molecules or "bions" have to do with localiza- 
tion [37|. A better understanding of the physical mech- 



anism behind localization in QCD could possibly shed 
some more light on the finite temperature chiral and de- 
confining transition. 
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